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We report on a quantitative investigation of the spin-dependent quasiparticle lifetimes and electron 
correlation effects in ferromagnetic hep Co(0001) by means of spin and angle-resolved photoemis- 
sion spectroscopy. The experimental spectra are compared in detail to state-of-the-art many-body 
calculations within the dynamical mean field theory and the three-body scattering approximation, 
including a full calculation of the one-step photoemission process. From this comparison we con- 
clude that although strong local many-body Coulomb interactions are of major importance for the 
qualitative description of correlation effects in Co, more sophisticated many-body calculations are 
needed in order to improve the quantitative agreement between theory and experiment, in particular 
concerning the linewidths. The quality of the overall agreement obtained for Co indicates that the 
effect of non-local correlations becomes weaker with increasing atomic number. 
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I. INTRODUCTION 



Over the last decades, the electronic structure and the 
dynamics of electronic states in solids have attracted a 
lot of attention. A better understanding of the binding 
energies and quasiparticle lifetimes in correlated systems 
like 3d ferromagnets- , transition metal oxides 2 -, 4f rare- 
earths 3 or high-T c superconductors 4 - has been achieved, 
together with important experimental progress in pho- 
toelectron and related spectroscopies^ On the theoret- 
ical side, the application of density functional theory 
(DFT) in the local density approximation 6 (LDA) has 
contributed with numerous calculations of single-particle 
E(k) band dispersions of solids, surfaces, and ultrathin 
films. 

In the case of magnetic transition metals, an exten- 
sive work in this respect has been carried out over the 
past thirty yearsi 7 - - — However, the question for the ex- 
act role of correlation effects in the bandstructure of 
these systems is still open and no general consensus has 
been reached due to several causes. The most impor- 
tant reasons for the disagreement between theory and 
experiments in these metals have been (i) the narrow- 
ing of the 3d bands due to correlation effects and (ii) 
the existence of new features in the photoemission data 
such as photoemission satellites in the core-level and 
valence band spectra which are not well described by 
single-particle approaches such as the generalized gra- 
dient approximation 16 (GGA) or the local spin-density 



approximation (LSDA) in DFT. 

Although one would expect a better agreement be- 
tween theory and experiment if the atomic number in- 
creases within the 3d series, the case of Co, which is the 
focus of the present work, is still an intriguing question. 
First attempts to compare experiment and theory indi- 
cated that for Co -but also for Fe- many-body effects can 
be considered smalliS, in contrast to Nii 17 i 18 ' 20 i 21 While 
Ni was considered an example of a strongly correlated 
system where the narrowing of the 3d bands is substan- 
tial, first results in Co indicated larger deviations than 
for Fo 10 ' 19 ' 22 where the energy positions with respect to 
the Fermi energy were about 10% smaller in the exper- 
iment than in LDA. 23 Moreover, the Fe photoemission 
peaks showed a large linear broadening (^60% of the 
binding energy)^ 3 - later on predicted to be accompanied 
by a strong loss in spectral weight 24 and most recently 
attributed to the strong non-locality of correlation effects 
in Fe^ 5 - On the other hand, Co has been the subject of 
several photoemission studies since the beginning of the 
80 a 10 ' 19 ' 26 " —, but only a few have been devoted to a de- 
tailed comparison between many-body calculations and 



experiments 



19.28.32 



Recently, by comparing experimental 



photoemission spectra of fee Co to theoretical calcula- 
tions within the three-body scattering approximation 33 
(3BS), it has been found that due to many-body ef- 
fects no sharp quasiparticle peaks exist for binding en- 
ergies larger than 2 eV in this system^ Interestingly, 
this effect has been theoretically explained by the exis- 
tence of strong correlation effects which could particu- 
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larly affect majority spin electrons 3 - 2 ., leading to a more 
pronounced renormalization of the majority spin quasi- 
particle spectral weight as compared to its minority coun- 
terpart. However, these conclusions were not supported 
by experiments involving spin resolution. On the con- 
trary, in a very recent photoemission experiment on a hep 
Co bulk single crystal with out-of-plane magnetization 34 
it has been shown that for sufficiently high photon en- 
ergies (488-654 eV), the spectral peak widths become 
narrower than previously observed in low-energy exper- 
iments and, as a result, quasiparticle bands of weak in- 
tensity can be observed even at binding energies larger 
than 3 eV. This effect has been attributed to the sig- 
nificant reduction of the final-state broadening contribu- 
tion to the measured linewidths at this photon energy 
range. In addition, a detailed analysis of the linewidths 
of the photoemission peaks revealed excellent agreement 
between experiment and calculations within the dynam- 
ical mean-field theory^ (DMFT). Although this conclu- 
sion was supported for both majority and minority spin 
electrons, once again, no spin resolution was provided 
in these experiments to decompose overlapping contribu- 
tions of opposite spin. 

In contrast to earlier studies, we utilize in the present 
paper spin resolution in order to give an experimental and 
theoretical corroboration to previous findings. The sys- 
tem under study is an ultrathin hexagonal close packed 
(hep) Co film of about 20 monoatomic layers (ML) thick- 
ness with an in-plane magnetizationj 29 ' 36 We compare 
the results from spin- and angle-resolved photoemission 
(SARPES) experiments to theoretical DMFT and 3BS 
many-body calculations, including a complete calculation 
within the one-step model (ISM) of photoemission. Our 
main conclusion is that although strong local many-body 
Coulomb interactions play a major role for the qualita- 
tive description of correlation effects in Co, more sophis- 
ticated many-body calculations including non-local inter- 
actions are needed in order to improve the quantitative 
agreement between theory and experiment, in particular 
concerning the linewidths. Although the results are sim- 
ilar to the ones previously reported for Fe^ the better 
agreement obtained in the case of Co indicates a reduc- 
tion of the non- locality of correlation effects with increas- 
ing atomic number. 



II. EXPERIMENT 

Experiments have been performed at room tempera- 
ture with a hemispherical SPECS Phoibos 150 electron 
energy analyzer and linearly polarized undulator radia- 
tion at the UE112-PGM1 beamline at BESSY II. For spin 
analysis, a Rice University Mott-type spin polarimeter 
has been operated at 26 kVj^Z The angular resolution 
of the equipment is better than 1° and the average en- 
ergy resolution about 100 meV. The Co(0001) surface 
has been prepared on W(110) by deposition of 20 ML 
Co from a high purity wire and post- annealing during 5 
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FIG. 1: (color online), (a) Geometry of the experiment, (b) 
Six-fold symmetry observed in the LEED pattern of the hep 
Co(0001) surface, (c) Sketch of the SBZ of the hep lattice. In 
(a) and (c), the magnetization orientation is also indicated. 

min. at ~ 400 °C. The tungsten substrate was cleaned 
before deposition by several cycles of annealing in oxy- 
gen (1 x 10~ 7 mbar) at 1200°C followed by flashing at 
2200°C. The base pressure during the experiments was 
1-2 x 10~ 10 mbar. Figure 1(a) shows an overview of 
the geometry of the experiment. The linear polariza- 
tion of the incident photon beam has been horizontal or 
vertical, and its angle of incidence with respect to the 
surface normal ^45°. Hence, in the reference frame of 
the sample, the light had more p- or full s-character, re- 
spectively. The structural quality of the Co films has 
been checked by low-energy electron diffraction (LEED). 
Fig. 1(b) shows the measured LEED pattern of the hep 
Co(0001) surface acquired with an incident electron beam 
of ^150 eV. Note that each diffraction spot corresponds 
to a r point of the surface Brillouin zone (SBZ), sketched 
in Fig. 1(c). The LEED pattern of Fig. 1(b) confirms 
that the annealing procedure which follows the epitaxial 
growth leads to a well-ordered surface of very high qual- 
ity. It is well-known that in hep Co(0001), the magnetiza- 
tion orientation lies out-of-plane in bulk single crystals 3 - 4 - 
and in-plane due to surface anisotropy effects in ultra- 
thin epitaxial films grown on W(110) i 29 ' 36 Therefore, in 
our experiments the sample was remanently magnetized 
in the film plane along the Co[1100] (or W[110]) easy 
axis^&^f^ 9 ., ag i nc ii ca ted in Fig. 1(a) and Fig. 1(c). Note 
that this direction is parallel to the T-M symmetry line 
of the SBZ. 



III. THEORY 

We start our considerations by a discussion of Pendry's 
formula for the photocurrent which defines the one-step 
model of photoemission 4 - : 

7 PES oc Im (ey.kiilG^AGf A f Ga |e/,k||> . (1) 

The expression can be derived from Fermi's golden 
rule for the transition probability per unit time. 41 
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Consequently, / PES denotes the elastic part of the 
photocurrent. Vertex renormalizations are neglected. 
This excludes inelastic energy losses and correspond- 
ing quantum-mechanical interference terms42r— Fur- 
thermore, the interaction of the outgoing photoelectron 
with the rest system is not taken into account. This 
"sudden approximation" is expected to be justified for 
not too small photon energies. We consider an energy-, 
angle- and spin-resolved photoemission experiment. The 
state of the photoelectron at the detector is written as 
|ey,k||), where k|| is the component of the wave vector 
parallel to the surface, and e/ is the kinetic energy of 
the photoelectron. The spin state of the photoelectron 
is implicitly included in |e/, k||) which is understood as a 
four-component Dirac spinor. The advanced Green func- 
tion G 2 in Eq. ([T]) characterizes the scattering proper- 
ties of the material at the final-state energy E 2 = ej- 
Via \^ f) = G^~|e/,k||) all multiple-scattering corrections 
are formally included. For an appropriate description 
of the photoemission process we must ensure the cor- 
rect asymptotic behaviour of ^/(r) beyond the crystal 
surface, i.e. a single outgoing plane wave characterized 
by ef and k||. Furthermore, the damping of the final 
state due to the imaginary part of the inner potential 
Voi(i?2) must be taken into account. We thus construct 
the final state within spin-polarized low-energy electron 
diffraction (SPLEED) theory considering a single plane 
wave |e/,k||) advancing onto the crystal surface. Us- 
ing the standard layer-KKR methoo— generalized for 
the relativistic cas o 44 i 45 , we first obtain the SPLEED 
state U*S? f (r) . The final state is then given as the time- 
reversed SPLEED state (U — —ia y K is the relativistic 
time inversion). Many-body effects are included phe- 
nomenologically in the SPLEED calculation, by using 
a parameterized, weakly energy-dependent and complex 
inner potential Vo(E 2 ) — V 0r (E2) + iV i(E 2 ). The imag- 
inary part Voi(E 2 ) was chosen constant in energy. The 
best agreement between theory and experiment was ob- 
tained for a value of 2.0 eV for all excitation energies, 
which means it was sufficient to take into account 20 
layers for the photocurrent calculation. For the quan- 
tum number 1 we used a maximum value l ma x=3 to ac- 
count for the contribution of d-f like transitions in the 
photoemission matrix elements^ This generalized inner 
potential takes into account inelastic corrections to the 
elastic photocurrent 4 ^ as well as the actual (real) inner 
potential, which serves as a reference energy inside the 
solid with respect to the vacuum level4£ Due to the finite 
imaginary part Vqi(E 2 ), the flux of elastically scattered 
electrons is permanently reduced, and thus the amplitude 
of the high-energy wave field VP / (r) can be neglected be- 
yond a finite distance from the surface. In order to take 
into account strong electronic correlation effects in the 
initial states we use the LSDA+DMFT approach realized 
in the framework of the fully relativistic Korringa-Kohn- 
Rostoker multiple scattering theory 48 (SPRKKR). The 
corresponding self-energy T, DMFT (E) is calculated fully 
self-consistently (e.g. in charge and self energy) via the 



DMFT approach^ Some of the results have been cross 
checked by non-self consistent LSDA+DMFT calcula- 
tions, implemented within the full potential linear muffin 
tin orbital (FP-LMTO) method^ This allows in particu- 
lar to use the full-potential option which might be impor- 
tant for more open structures. As a DMFT-solver the rel- 
ativistic version of the so-called Spin-Polarized T-Matrix 
Plus Fluctuation Exchange (SPTF) approximatio n 50 ! 5 1 
was used. In the present study, all the calculations were 
performed with a bulk-like self-energy, i.e. not layer- 
dependent. However, corresponding Bloch spectral func- 
tions were cross-checked with the FP-LMTO results us- 
ing layer-dependent self-energy and only weak quantita- 
tive differences were found. 

In addition to these approaches we accounted for 
correlation effects within the three-body scattering 
approximation^ Within the 3BS-approach the self- 
energy is calculated using a configuration interaction-like 
expansion. In particular three-particle configurations like 
one hole plus one electron-hole pair are explicitly taken 
into account within 3BS-based calculations. The corre- 
sponding output can be directly related to the photoe- 
mission process and allows for a detailed analysis of var- 
ious contributions to the self-energy (e.g., electron- hole 
lifetime) . 

IV. RESULTS AND DISCUSSION 

In this section, we discuss the effects of electron corre- 
lations in the spin-dependent electronic structure of hep 
Co. In subsection A, the relativistic notation used to 
label the symmetries of the observed states will be in- 
troduced. In subsection B, a detailed analysis of the 
different transitions observed in the experimental and 
theoretical SARPES spectra will be given, focusing on 
the quality of the agreement between theory and exper- 
iment. Finally, in subsections C and D, further discus- 
sion on the many-body aspects of the Co bandstructure 
will be presented, together with a quantitative analysis 
of the experimental and theoretical linewidths of the Co 
photoemission peaks. 



A. Selection rules and hep notation 

In contrast to previous studies of hep Co^ where the 
non-relativistic notation for fee crystals was used to la- 
bel the symmetry of the observed states, from now on we 
will follow the relativistic notation given by Benbow 52 
for hep crystals. This notation is much simpler since 
due to the lower symmetry it only involves initial and 
final states of A7, As and Ag symmetries. This means 
that in the geometry used in the experiment (Fig. 1(a)) 
and due to selection rules, only transitions of A7 and Ag 
symmetries are suppressed for s-polarization when mov- 
ing from r to A and from A to T in the next Brillouin 
zone, respectively. Hence, only transitions of Ag sym- 
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FIG. 2: (color online), (a) Experimental spin- integrated photoemission spectra of hep Co(0001) measured with p-polarization 
in normal-emission along the TA direction of the bulk Brillouin zone at different excitation energies (steps of 2 eV). The curves 
are labeled by selected photon energies and by the wave vectors in units of rA=0.77A -1 (r point at x=0 and A point at 
x=l). (b) Corresponding calculation obtained by the LSDA+DMFT+1SM method for hep Co with an in-plane magnetization 
along the [1100] direction, (c) and (d): Spin-integrated band structure calculations for two configurations of the magnetization, 
out-of-plane (along [0001] direction, green (light) lines) and in-plane (along [1100] direction, black (dark) lines). In (c), LSDA 
calculations are shown, and in (d) fully relativistic valence band states obtained self-consistently in the SPRKKR formalism 
by suppressing the imaginary part of the self-energy IhiE^mft- 



metry can be observed with s-polarized photons. The 
effect of matrix-elements, which is also influenced by the 
symmetry of the hep crystal structure^ 3 -, is similar to 
the interference effect or Brillouin zone-selection rule re- 
ported in Ref. l54l which suppresses the intensities of 
certain bands in a particular Brillouin zone but increases 
them in the next-higher one. The overall effect leads to 
a bandstructure which is equivalent to the one of an fee 
crystal if all the hep bands are folded back about the 
A points This also implies that states of A 7 and A§ 
symmetry should be interchanged when crossing the A 
point while states of Ag symmetry remain unaffected. 
To avoid redundancy, from now on double subscripts will 
be adopted, unless specified. For better clarity, the sub- 
script notation of almost degenerated states will follow 
their order of appearance from lower to higher binding 
energy (BE). 



B. Spin and angle- resolved photoemission: a 
comparison theory-experiment 

Figures 2(a) and 2(b) show a comparison between 
the spin-integrated experimental spectra of hep Co(0001) 
and the corresponding theoretical LSDA+DMFT+1SM 



calculations for p-polarization along the TA direction of 
the bulk Brillouin zone. The k values were calculated 
from the measured photon energies, ranging from 22 to 
80 eV and using an inner potential Vb=14.8 eV. The best 
correspondence between the BE positions of the exper- 
imental and theoretical peaks was found for values of 
the averaged on-site Coulomb and exchange interactions 
of U=2.5 eV and J=0.9 eV, in agreement with previous 
theoretical and experimental studiesi 32 ' 35 Therefore, un- 
less specified, these have been the parameters of choice 
in all the calculations presented here. 

Figures 2(c) and 2(d) show LSDA and fully-relativistic 
LSDA+DMFT calculations of the spin-integrated band 
structures for different directions of the sample magneti- 
zation, the out-of-plane [0001] and in-plane [1100] direc- 
tions (green (light) and black (dark) lines, respectively). 
It can be observed how the different bands shift between 
Fig. 2(c) and Fig. 2(d) towards lower binding energies 
due to electron correlation effects for both magnetization 
orientations. 

As already mentioned above, the Co(0001) films are 
magnetized in-plane along the [1100] direction. All 
LSDA+DMFT+1SM calculations presented here have 
been performed using in-plane magnetization direction. 
However, we show in addition out-of-plane calculations 
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in order to emphasize how due to spin-orbit coupling 
effects, changes in the orientation of the magnetization 
correspond to changes in the electronic structure of hep 
Co. This purely relativistic effec t 45 ' 55 i 56 may gain im- 
portance if e.g., experimental and theoretical spectra of 
hep Co bulk single crystals 2 -^ and overlayer a 29 ' 32 need 
to be compared. To a first approximation, the magne- 
tocrystalline anisotropy energy is given by the average 
energy difference between the BE positions of the bands 
in Fig. 2(c) and Fig. 2(d) with in-plane and out-of- 
plane magnetization directions, respectively. The overall 
effect of the magnetic anisotropy on the bandstructure 
results in slight changes in the dispersion behaviour of 
some bands, additional BE shifts and most importantly, 
a considerable reduction of the spin-orbit splitting with 
an in-plane magnetization configuration. In this case, 
near the Fermi level and below 2 eV, most of the bands 
are degenerated while in the out-of-plane configuration, 
the spin-orbit splitting is ~0.1 eV in average, a value of 
the order of the experimental energy resolution. There- 
fore, it should be emphasized that relatively narrower 
peaks would be expected in photoemission experiments 
with an in-plane configuration of the magnetization, such 
as in Co overlayers as opposed to hep Co bulk single crys- 
tals. 

Let us now focus on the origin of the different transi- 
tions appearing in the theoretical and experimental spec- 
tra of Fig. 2 and Fig. 3, from lower to higher binding 
energies. In Fig. 3, similar data to the ones presented 
in Fig. 2(a) and Fig. 2(b) are shown, but now at a few 
selected photon energies and spin resolved. Figures 3(a) 
and 3(b) show experimental SARPES spectra for p- and 
s-polarized photons, respectively, while Fig. 3(c) and Fig. 
3(d) show the corresponding LSDA+DMFT+1SM calcu- 
lations. All the peaks appearing near the Fermi level in 
the full photon energy range and for a BE lower than 
Eb ^0.5 eV, can be attributed to the two majority and 
minority spin components of a Tamm-like surface feature, 
in agreement with previous experimental observations. 31 
Closest to the Fermi level, the minority spin component 
of this state appears at Eb ^0.05 eV, and it can be iden- 
tified as a pure surface state of Ag symmetry, since it lies 
within the gap appearing near the Fermi level and its 
intensity is nearly the same for both p- and s-polarized 
photons (see Fig. 3). Only near the first Y point and at 
x~0.04rA (h^=26 eV), this peak is slightly influenced 
by the minority spin bulk-like band of Ag 7 symmetry 
crossing the Fermi level. Thus, only at this photon en- 
ergy it acquires a slight bulk-like character and it can be 
identified as a surface resonance feature. 

A similar argumentation can be used for the major- 
ity spin component of the abovementioned surface-like 
Tamm state, which can be identified as a surface res- 
onance feature of Ag 7 symmetry (or equivalently, Ag 8 
symmetry in the next Brillouin zone), since it is located 
at the border of the gap and almost degenerated with 
the majority bulk bands appearing around the same BE. 
This state, which appears at Eb ~0.4 eV, can be ob- 
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FIG. 3: (color online). Spin-resolved spectra of hep Co(0001) 
for different polarizations, (a), (b) Experiment (upwards 
(black (dark)) triangles: majority spin states; downwards 
(red (light)) triangles: minority spin states). (c), (d): 
LSDA+DMFT+1SM theory (black (dark) and red (light) 
lines for majority and minority spin electrons, respectively). 
(a),(c) for p- and (b), (d) for s-polarization. Additional up- 
wards (black) [downwards (red)] triangles on top of the spec- 
tra indicate majority [minority] bulk states, while vertical 
black (red) bars majority (minority) spin surface-related fea- 
tures. The curves are labeled by the photon energies (right) 
and the wavevectors (left) for both the hep (rA=0.77A _1 ) 
and fee (FL=1.54A -1 ) lattices for comparison. 



served as a non-dispersing feature in the full photon en- 
ergy range (see Fig. 2(a), Fig. 2(b) and Fig. 3). How- 
ever, when moving in the experiment from Y to A in a 
photon energy range between 26 and 38 eV (0.04< x < 
0.58) (Fig. 2(a)), its intensity is much weaker as com- 
pared to the theoretical calculations (Fig. 2(b)). This 
can also be noticed in the results of Fig. 3 near the Y 
point (x~0.08 TA, hv—26 eV) for both p and s polariza- 
tions. In this case it can be observed as a majority spin 
shoulder close to the Fermi level, but in the theoretical 
spectra as a sharp peak. Since in the calculations its in- 
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tensity does not seem to be strongly influenced by matrix- 
element effects and in addition it exhibits a significantly 
narrower linewidth, we attribute these strong differences 
to a theoretical underestimation of the multiple scatter- 
ing events occurring between the surface and the bulk 
electron wavefunctions. Moreover, the intensity of this 
peak is strongly enhanced at higher energies in the ex- 
periment (see Fig. 2(a), Fig. 3(a) and Fig. 3(b)), partic- 
ularly when moving from around the A point (x— 0.98rA, 
lw=48 cV) to the second T point (x-0.0irA, hv=72 eV) 
and beyond. This increase in intensity can be attributed, 
on the one hand, to a slightly reduced surface sensitiv- 
ity at higher photon energies, and on the other hand, to 
the contribution of the nearly-flat majority spin bulk-like 
band of Ag 8 symmetry appearing around the same BE 
for photon energies hi/ >48 eV. Note once again, that due 
to the symmetry arguments mentioned above, for hi/ <48 
eV this bulk band has Ag 7 symmetry and is located at 
a BE of about 0.7 eV near the first T point (hi/=26 eV) 
(Fig. 3). 

All the other bands that are left to describe are bulk- 
like and appear at a BE higher than —0.5 eV for all 
the measured k points. Firstly, we observe a minority 
spin band of Ag 7 symmetry located in the experiment at 
E B -0.8 eV in the first T point (see Fig. 2(a), Fig. 3(a) 
and Fig. 3(b)). Through the polarization dependence 
of the experimental spin-resolved data of Fig. 3(a) and 
Fig. 3(b), it is possible to see that this band comprises 
several components which are all degenerated at the first 
r point. These components slightly disperse when the 
photon energy is increased, giving rise to shoulders which 
can be detected through polarization-dependent analysis. 
Most prominently, at x=0.67rA (h^=40 eV) a minority 
spin shoulder at a BE of about 0.75 eV can be clearly 
distinguished in the lower BE energy side of the peak 
appearing at Eb —0.95 eV for p-polarization data (Fig. 
3(a)). Detailed analysis reveals that the BE position of 
this shoulder exactly coincides with the minority spin 
peak of Ag symmetry appearing in the corresponding 
spectrum for s-polarization (Fig. 3(b)). Further infor- 
mation on the symmetry of these states can be obtained 
from a comparison to the spin-resolved theoretical spec- 
tra of Fig. 3(c) and Fig. 3(d). In particular, it can be 
demonstrated that the peak appearing at 0.95 eV must 
be assigned to a state of A 7 symmetry, and the shoulder 
to a state of A 7 9 symmetry. A similar effect occurs at 
higher photon energies in the next Brillouin zone e.g., at 
x=0.62rA (hi/=56 eV), where the peak at E B -1.2 eV 
and the shoulder at Eb —0.8 eV can be assigned to bands 
of Ag and Ag 9 symmetries, respectively. In the next T 
point and at x=0.08rA (hi/=70 eV), the peaks show the 
same behaviour, with the difference that a weak shoul- 
der appears at Eb —1.5 eV, exactly in the higher BE 
side of the bulk-like Ag 9 state located at Eb —0.8 eV. 
Only by means of a detailed comparison to the theoreti- 
cal spectra of Fig. 3(b), we can attribute this feature to 
a Ag g surface resonance. The fact that it was not visible 



at lower photon energies is possibly related to the dis- 
persion of the bulk-like minority spin peaks we have just 
discussed. Besides, this may also be related to the very 
short lifetime of surface resonance features at higher BE 
caused by the strong increase of IuiSdmft, which also 
leads to a very large broadening and makes it impossible 
to properly distinguish them in the experiment. 

At last, we are only left with the majority spin bulk-like 
states at higher BE. These peaks are strongly broadened 
and more difficult to distinguish in the experiment, par- 
ticularly in the spectra of Fig. 2(a) due to the scaling 
of the intensities. However, with the use of spin resolu- 
tion (see Fig. 3(a) and Fig. 3(b)), it can be seen that 
they are mostly suppressed with s-polarization at bind- 
ing energies around 2 eV. These states become clearly 
visible in a BE range from about 1.4 to 1.9 eV for p- 
polarization data in Fig. 3(a) at x=0.67rA (hi/=40 eV) 
and x=0.62rA (hz/=56 eV), and for s-polarization data 
in Fig. 3(b) near the first T point at x=0.04rA (hi/=26 
cV). In the calculations (Fig. 3(c) and Fig. 3(d)), on 
the other hand, the BE of these states seems to be well 
reproduced, but they are more pronounced in intensity. 
From the polarization dependent analysis of the exper- 
imental data and by comparison to the theoretical cal- 
culations, the assignment of the symmetries in this case 
can be derived from selection rules in a straightforward 
manner. Near the first T point eV), where all the 

states of this kind are degenerated, we can observe them 
with both polarizations, and thus the dominant symme- 
try should be A|. At x=0.67rA (hz/=40 eV), however, 
the dominant symmetry should be A 7 , since the peaks 
are mostly suppressed with s-polarization. Moving into 
the next Brillouin zone, the same argument applies and 
therefore, the dominant symmetries are Ag and Ag 9 for 
x=0.62rA (hz/=56 eV) and 0.08rA (hz/=70 eV), respec- 
tively. 

Let us now discuss the quality of the agreement be- 
tween the theoretical and experimental spectra shown in 
Fig. 2 and Fig. 3. The first observation is that we find 
rather good agreement for most of the BE positions of 
the peaks appearing in the spectra. The second observa- 
tion, concerning the linewidths, is that the larger broad- 
ening of the experimental peaks indicates that the theory 
underestimates the scattering rates, similar to what has 
been found in Fe^ The third observation, in contrast 
to recent experimental conclusions^, is that no quasi- 
particle bands derived from 3d states appear at binding 
energies higher than —2 eV. At around 2 eV, this effect 
is even more pronounced in the experiment than in the 
theory due to the increased linewidth broadening men- 
tioned above. This finding is in agreement with previous 
experimental and theoretical analysis in the framework 
of the 3BS theory^, in which the quenching of the major- 
ity spin channel quasiparticle excitations was identified 
as the main mechanism for this intriguing effect. Our 
results corroborate this conclusion in view of the spin- 
resolved experimental and theoretical spectra presented 
in Fig. 3. Only peaks with strong sp character are ob- 
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served at higher BE energies than 2 eV, such as the mi- 
nority spin sp-band of Ag symmetry appearing at x^0.08 
TA (h^=70eV) for p-polarization in Fig. 3(a) and in Fig. 
3(c) at a BE of about 4.3 eV. 



C. Many-body aspects of the bandstructure 

Comparing the experimental results from spin- 
integrated and spin-resolved ARPES measurements with 
LSDA+DMFT+1SM results, good agreement is obtained 
for most of the peak positions. This is also shown in Fig. 
4, where the BE positions of the different experimental 
peaks we have just described arc compared to the calcu- 
lated bulk-like spin-dependent spectral functions A(k,E) 
obtained in the framework of the LDSA+DMFT (Fig. 
4(a) and Fig. 4(b)), LSDA+3BS (Fig. 4(c) and Fig. 
4(d)) and LSDA (Fig. 4(e) and Fig. 4(f)) schemes. 
Majority spin states are shown in Fig. 4(a), Fig. 4(c) 
and Fig. 4(e) and minority spin states in Fig. 4(b), 
Fig. 4(d) and Fig. 4(f). The symmetry labels of the 
states are also given. In these calculations, the experi- 
mental energy resolution is considered by a broadening of 




r a r a r a 

FIG. 4: (color online). Quasiparticle spectral functions 
of hep Co and the photoemission peak positions obtained 
from the spin-resolved measurements of Fig. [3] for differ- 
ent polarizations (O for p- and O for s- polarization). The 
results are obtained by the LSDA+DMFT ((a),(b)), the 
LSDA+3BS ((c), (d)) and the LSDA ((e),(f)) methods for ma- 
jority ((a),(c),(e)) and minority spin ((b),(d),(f)) electronic 
states, respectively. The symmetries of the different states 
are indicated. 



^100 meV. In a first approximation, Figs. 4(a)-4(d) give 
a comparison between the spin-dependent experimental 
ReE and theoretical ReT,DMFT and ReT, 3 BS- However, 
it should be noted that the spectral functions A(k,E) 
contain more bands than the ones observed in the ex- 
periment, since Brillouin zone selection effects are not 
considered in this case. This means that in these cal- 
culations matrix-element effects, surface effects, and the 
actual steps of the photoemission process (i.e., propaga- 
tion of the photoelectron to the crystal surface and its es- 
cape into the vacuum) are not included. The experimen- 
tal peak positions were extracted from the experimental 
SARPES data of Fig. 3(a) and Fig. 3(b) by means of a 
fitting procedure which will be described in more detail 
below. Before proceeding, it should be noted that from 
Figs. 4(a)-4(d), it is clear that the theoretical results de- 
rived from both LSDA+DMFT and LSDA+3BS methods 
are well in agreement with each other for identical values 
of U and J. 

Let us focus first on majority spin states. From a first 
examination across the complete Brillouin zone, one can 
conclude that the LSDA+DMFT and LSDA+3BS spec- 
tral functions reproduce the experimental data better 
than standard LSDA calculations. This is most apparent 
at lower BE and for states near the A point, where the 
LSDA bands are shifted to higher BE as compared to 
the experimental points assigned to majority spin states 
of Ag 8 symmetry, as described above. Thus, in order to 
find best agreement with the experiment, we would need 
to shift by about 0.5 eV the corresponding majority spin 
LSDA bands towards lower BE. Moreover, note the re- 
markable intensity of the LSDA majority spin bands ap- 
pearing in a BE range from 2.5 to 3 eV in the whole Bril- 
louin zone. As discussed before, these bands are strongly 
broadened in the experiment and appear at BE energies 
of around 1.7-1.9 eV. From the polarization-dependent 
analysis carried out before, we also know that at the 
different k-points of the bandstructure, the symmetries 
of these states [Ag(As.g) and Ay(Ag)] coincide with the 
symmetries of the LSDA bands appearing at Eb >2 eV, 
which are not forbidden by selection rules. This means 
that in Fig. 4(e), the agreement between the peak posi- 
tions from the experimental peaks around 1.9 eV with the 
LSDA bands of Ag 7 appearing at binding energies below 
2 eV is just accidental, since this agreement does not 
hold based on simple symmetry arguments. Therefore, 
also in this case a shift of these states by about 0.7-1 eV 
becomes necessary to properly explain the bandwidth re- 
duction observed in the experiment. These two facts are 
already an indication that a narrowing of the 3d bands 
due to correlation effects becomes necessary to properly 
describe the experimental data. This is clearly demon- 
strated in Fig. 4(a) and Fig. 4(c), where a relatively 
good agreement is obtained for the BE positions of the 
experimental majority spin peaks for both p and s po- 
larizations. In this reduction of the bandwidth of 
majority spin states by ~37% is necessary, in agreement 
with the results presented in Fig. 2(c) and Fig. 2(d). 
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Regarding the minority spin states, it is clear that cor- 
relation effects are weaker than for majority spin. Note 
that close to the A point, the Ag minority spin compo- 
nent of the surface-like Tamm state discussed above is 
also visible. This state is not reproduced in the spectral 
function calculations because, as already explained, sur- 
face effects are not included in these calculations. In fact, 
the different minority spin bands in the LSDA+DMFT 
(Fig. 4(b)) and LSDA+3BS (Fig. 4(d)) calculations 
shift by about 0.1 eV respect to the LSDA ones (Fig. 
4(e)). This means that the BE positions of the differ- 
ent bands in the three different theoretical approaches 
show relatively good agreement with the experimental 
data. This is something expected a priori, since from the 
theory we know that correlation effects should be more 
pronounced for majority spin electrons than for minority 
spin ones, which is a common feature for 3d transition 
metals. This can be easily explained by considering the 
creation of electron-hole pairs in the minority spin chan- 
nel: since there are less minority spin electrons and thus 
more empty states of this kind, for any processes involv- 
ing electron-hole pair creation, the pair is more likely 
to appear in the minority spin bands. This means that 
any scattering process involving majority spin electrons 
mostly leads to the creation of minority spin electron- 
hole pairs, with an effective interaction in the complete 
process proportional to U. On the other hand, scattering 
processes involving minority spin electrons also lead to 
the creation of minority spin electron- hole pairs, but the 
effective interaction for parallel-spin pairs is proportional 
to U-J<U. Therefore, correlation effects are stronger for 
majority spin electrons and no big differences are ob- 
served between the results presented in Fig. 4(b), Fig. 
4(d) and Fig. 4(f). 

The next point that we should discuss is the 
agreement between the experimental and theoretical 
LSDA+DMFT+1SM spectra as a function of the on-site 
Coulomb interaction U. As already mentioned, the best 
agreement between theory and experiment was found 
for U=2.5 eV and J=0.9 eV. The value of J adopted 
here can be considered a good general assumption for 
all 3d elements, since J is only a weakly screened pa- 
rameter. This means that no differences are observed in 
the LSDA+DMFT+1SM spectra for different J values. 
Therefore, the only parameter directly linked to the ef- 
fect of electronic correlations which might lead to changes 
in the calculated spectra is U and we will proceed with 
discussing how sensitive the theory is to this parameter 
in view of the experimental results. 

In Fig. 5, we compare the spin-integrated exper- 
imental spectra measured for p-polarization at the T 
point (x=0.0irA, hi/=72 eV) [Fig. 5(a)] and near the 
A point (x=0.96rA, hz/=48 eV) [Fig. 5(c)] to the 
LSDA+DMFT+1SM and LSDA calculations broadened 
by the experimental energy resolution. In this case, the 
LSDA calculations also include surface effects for a bet- 
ter comparison. Among the two spectra available in the 
data near the F point, the spectrum at Iw=72 eV has 
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FIG. 5: (color online). Comparison between experimental 
and theoretical spectra at the T [(a) and (b), hv=72 eV] 
and A [(c) and (d), hz/=48 eV] points for p-polarization. (a) 
Spin-integrated experimental spectra [thick (black) symbols], 
single-particle LSDA-based calculation including surface ef- 
fects [thick dotted (green) line] and LSDA+DMFT+1SM 
spectra [thick solid (green) line]. (b) Spin-dependent 
LSDA+DMFT+1SM calculations for different U values. 
Black (dark) [red (light)] color denote majority [minority] spin 
and green (lighter) color the corresponding spin-integrated 
spectrum. Thin dotted, dashed, thick and thin solid lines for 
U=1.5, 2, 2.5 and 3 eV, respectively. In (c) and (d), analogous 
data as in (a) and (b). Symmetry labels in the two consecu- 
tive Brillouin zones of the hep lattice are also indicated. 



been selected as a suitable choice because of the better 
agreement at higher photon energies between the exper- 
imental and theoretical intensities of the majority spin 
Ag g surface resonance located near the Fermi level. The 
spin character and symmetry of the different experimen- 
tal states are also given. In Fig. 5(b) and Fig. 5(d), 
spin-dependent LSDA+DMFT+1SM calculations for U 
ranging from 1.5 to 3 eV are also shown. 

In general, the Co LSDA+DMFT+1SM calculations 
lead to an important improvement as compared to the 
LSDA calculations. In both Fig. 5(a) and Fig. 5(c), 
the BE positions of the theoretical peaks agree well 
with the experimental ones if correlation effects are in- 
cluded. Regarding majority spin states, in Fig. 5(a) a 
good agreement is achieved for the BE position of the 
Ag g surface resonance near the Fermi level for both the 
LSDA+DMFT+1SM and LSDA calculations. The bulk 
component of this peak, which appears at Eb —0.78 
eV in the LSDA calculations, is shifted due to corre- 
lation effects in the LSDA+DMFT+1SM calculation to 
E B -0.47 eV (see Fig. 5(b), U=2.5 eV). In the experi- 
ment these peaks are almost degenerated at that BE due 
to the experimental energy resolution. A similar situa- 
tion can be found in Fig. 5(c), where good agreement 
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is obtained in the BE position of the Ag 7 majority spin 
component of the Tamm-like surface state. 

Let us focus now on the most prominent bulk-like mi- 
nority spin peaks appearing in the theoretical and exper- 
imental spectra. The Ag 9 minority spin bulk-like state, 
which in the LSDA calculation is located at a E B —0.8 
eV, is shifted to E B —0.66 eV in the experiment (Fig. 
5(a)) and to E B -0.57 eV in the LSDA+DMFT+1SM 
calculation (Fig. 5(b)) due to correlation effects. The 
Ay g minority spin bulk-like state, on the other hand, is 
shifted to Eg —0.79 eV in the experiment (Fig. 5(c)) 
and to E B -0.63 eV in the LSDA+DMFT+1SM cal- 
culation (Fig. 5(d)) with respect to the LSDA value of 
E B -0.96 eV. At last, the A^ peak is shifted to E B -1.20 
eV in the experiment (Fig. 5(c)) and to E B —1.22 
eV in the LSDA+DMFT+1SM calculation (Fig. 5(d)) 
compared to E B —1.38 eV in the LSDA calculation. 
This picture is consistent with an averaged experimen- 
tal mass enhancement m* /mo wl.26, which should be 
compared to the theoretical value of m* /mo «1.29 for 
U=2.5 eV. A similar analysis can be carried out for the 
LSDA+DMFT+1SM spectra shown in Fig. 5(b) and Fig. 
5(d) for U values around 2.5 eV. In particular, we obtain 
m* /mo —1.22 for U=2 eV in average, and m* /mo «1.55 
for U=3 eV. Although the deviations between the experi- 
mental and theoretical mass enhancement factors for dif- 
ferent U values are not large, from this evaluation U=2.5 
eV is closest to the experiment. However, a precise de- 
termination of U requires a more sophisticated approach 
based on a full analysis of this type in the whole Brillouin 
zone. This was done in order to decide the best choice of 
U in the present calculations: the complete set of exper- 
imental and LSDA+DMFT+ ISM spectra were directly 
compared in terms of BE position, intensity and width 
of the peaks over the whole Brillouin zone for different 
values of U, and better overall agreement was found for 
U=2.5 eV. Another criterion which was also combined 
to the one just mentioned is related to the spin char- 
acter of the different bands. An example showing this 
can be noticed in e.g., Fig. 5(a) and Fig. 5(b), where 
the LSDA+DMFT+1SM spin-resolved spectra for U=2 
eV exhibit a majority spin peak at E B —0.6 eV which 
in turn, is at the same BE as the experimental minority 
spin peak of Ag 9 symmetry. 

Furthermore, it can be seen in Fig. 5(b) and Fig. 5(d) 
that increasing the value of U does not shift the peaks sig- 
nificantly towards the Fermi level. This can be attributed 
to the following reasons: (i) Most of the intensity in the 
spectra for BE<0.5 eV is due to surface-related features, 
which are not strongly sensitive to changes in the U pa- 
rameter, (ii) the shift of pure bulk-like states only occurs 
in a narrow energy interval from 0.5 to 1.8 eV, other- 
wise their sensitivity to U is reduced through coupling 
to surface-related states near the border of the gap, (iii) 
there are no strong changes in T,omft in the region very 
close to the Fermi level. Because of these reasons, in Co 
the peaks exhibit smaller shifts with increasing U than 



in, e.g. Fe, because they are closer to the Fermi level, 
and most of the bulk bands are located in a narrow en- 
ergy interval. Besides, the coupling between surface and 
bulk states also plays a role, since in Co there is a gap 
around the T point. This is also related to the fact that 
due to the hexagonal structure the bands are relatively 
flat, leading to real surface states which are less sensitive 
to changes in U and a weak bulk-surface coupling. 



D. Spin-dependent quasiparticle lifetimes 

The last issue to be discussed is the quality of the 
agreement between the experimental and theoretical 
linewidths of the Co photoemission peaks. As it was 
pointed out above, it is clear that the Co calculations 
presented here underestimate the scattering rates, simi- 
lar to what we found in Fe.— In Fig. 5, it can also be 
noticed that the width of the peaks does not increase with 
U. In the narrow BE range we are focusing on, increasing 
U shifts the peaks to lower BE and also leads to nearly- 
constant values of the energy-dependent IibSdmft- Be- 
sides, since the effect of U on the peak positions is rel- 
atively small, a change in the linewidth can hardly be 
appreciated. This holds for both bulk-like and surface- 
related features. 

We emphasize that the additional experimental broad- 
ening cannot be caused by final state effects^, since (i) 
their contribution to the linewidth of surface states can 
be neglected because no dispersion normal to the sur- 
face exists in this case, (ii) they can be neglected when 
compared to the experimental resolution near the criti- 
cal points (e.g. at the T and A points in Fig. 5), (iii) 
they are completely considered in the ISM calculations, 
(iv) they are expected to be small for Co(0001) because 
due to the hexagonal structure the bands do not show 
strong initial state dispersion. Because of these reasons, 
we exclude final state effects as possible mechanism for 
the extra broadening observed in the experiment. 

Nevertheless, in the following we would like to obtain 
more quantitative information on the discrepancy be- 
tween experiment and theory we have observed concern- 
ing the linewidths. As mentioned before, a fitting proce- 
dure was used to extract the experimental peak positions 
from SARPES experiments shown in Fig. 3(a) and Fig. 
3(b). This evaluation is also useful for an estimation by 
how much the theoretical linewidths are underestimated 
compared to the experimental results. The use of spin- 
resolution is remarkably important in this case because 
from the spin-dependent behaviour of the self-energy in 
the theoretical calculations one would expect the exper- 
imental linewidths of initial states to be spin-dependent 
as well. 

The experimental spin-resolved spectra were fitted at 
various k values by a sum of Lorentzians plus a back- 
ground. The fitting procedure works as follows: for a 
particular k-point, a SARPES spectrum containing N 
peaks is fitted by a function involving a convolution of 
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the form: 




f(E,T) 'Y^ M i- A ki (Ei,Ui) + B k {E)j ® G{hv) 



(2) 

where f(E,T) is the Fermi function and Ei, uji and 
the matrix-elements Mi are fitting parameters corre- 
sponding to the BE, width and intensity of the quasi- 
particle peaks for different polarizations, respectively. 
The spectral function Aki(Ei,uji) is approximated by 
Lorentzian functions and Bk (E) assumed to be a Shirley- 
like background^. Here the subindex k denotes the kj^ 
points. The full width at half maximum (FWHM) of the 
Gaussian slit function G(hv) corresponds to the total en- 
ergy resolution A.E to t of the experiment, which is photon 
energy- dependent, i.e.: 

AE tot = v/(A£ Q ) 2 + (AE hv )* (3) 

where AE a and AEh u are the energy resolutions of the 
electron energy analyzer and beamline monochromator, 
respectively. The Shirley-like background Bk(E), which 
accounts for the effect of the inelastic scattering of elec- 
trons, is calculated as: 

B k (E) = C-(j^ dE -J2 M ? ■ Akt{Ei,Ui)j (4) 

where C is a fitting constant and the integration is done 
over the complete spectrum with E as a running param- 
eter over the kinetic energy axis. 

Figure 6 shows few selected results of the fits obtained 
by this procedure for p (Fig. 6(a) and Fig. 6(b)) and 
s (Fig. 6(c) and Fig. 6(d)) polarized photons, and for 
both majority (Fig. 6(a) and Fig. 6(c)) and minority 
(Fig. 6(b) and Fig. 6(d)) spin states. We should keep in 
mind that such an evaluation is problematic due to the 
strong energy- dependent behaviour of ImS which leads 
to the formation of asymmetric Lorentzians. The appear- 
ance of such an asymmetry can be associated with the 
damping of quasiparticle excitations and the correspond- 
ing increase of the incoherent part of spectral function 
with increasing BE. However, at energies close enough 
to the Fermi level, we would expect this asymmetry to 
be strongly reduced since most of the peaks are due to 
coherent excitations which correspond to quasiparticlcs 
with well-defined energy and momentum. 

For the fitting procedure we have adopted an experi- 
mental approach: almost degenerated bands which can- 
not be resolved due to the experimental energy resolu- 
tion are fitted with a single peak. An exception to this 
point occurs when some bands can be suppressed due 
to matrix elements effects. In this context, polarization- 
dependent measurements are suited to accurately obtain 
the linewidth of specific bands that are nearly degener- 
ated and cannot be resolved independently under cer- 
tain light polarization conditions. This is the case of 
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FIG. 6: (color online). Fits to spin-resolved valence band 
spectra at x=0.62rA (hv=56 eV) of the hep Co(0001) bulk 
Brillouin zone, (a), (b) for p- and (c), (d) for s-polarization. 
Black (dark) [red (light)] lines are the results of the fits to ma- 
jority [minority] spin spectra, which are indicated by upwards 
[downwards] triangles and arrows, respectively. Thin lines 
correspond to the different Lorentzian functions and back- 
grounds used in the fits. Thick lines indicate the final fit of 
each spectrum obtained by using Eq. ([2]). Symmetry labels, 
surface resonance (SR) and surface state (SS) features are 
indicated. 



e.g., the majority spin states of Ag 9 symmetry located 
at Eb ~1.4 eV for p-polarization data (Fig. 6(a)), which 
form a shoulder in the corresponding s-polarization spec- 
tra around the same BE (Fig. 6(b)). The same holds 
for the Ag 9 minority spin states located at Eb ~0.7 eV 
for s-polarization data (Fig. 6(d)), which in turn form 
a shoulder in the corresponding p-polarization spectra 
(Fig. 6(b)). This allows us to fix the width and BE po- 
sitions obtained from the fits of individual bands for a 
given polarization in order to fit the shoulders of almost 
degenerated bands which are observed with the opposite 
polarization. 

Concerning the lincwidths, the result of this analysis 
is presented in Fig. 7, where we have extracted the spin- 
dependent experimental lmTi exp of Co as a function of 
BE and k and compared it to the theoretical one. The 
results for majority (minority) spin electrons are shown 
in Fig. 7(a) (Fig. 7(b)), together with the correspond- 
ing LSDA+DMFT (U=2.5 and 3 eV) and LSDA+3BS 
calculations (U=2.5 eV). It is clear that no quantitative 
agreement between the theoretical calculations and the 
experimental data is observed. The experimental data 
points have been corrected by all the non-electronic con- 
tributions to the broadening with the exception of im- 
purity scattering (Im£j mp = Ti mp /2, where r=FWHM) 
since it can only be determined once all the other contri- 
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FIG. 7: (color online). Comparison between the experimental 
(symbols) and theoretical (lines) imaginary part of the self- 
energy of hep Co(0001) for (a) majority (black (dark) color) 
and (b) minority (red (light) color) spin electrons. The exper- 
imental data points only contain the electronic and impurity 
scattering contributions to the linewidth. On the right hand 
side, labels to identify the different k-points are also given for 
both hep and fee lattices. Full and open symbols are for p- 
and s-polarization, respectively. The theoretical calculations 
correspond to ImSoMFT for U=2.5 eV (thick dotted lines) 
and U=3 eV (thin dotted lines) and to lmS^BS f° r U=2.5 eV 
(solid lines). 



butions have been corrected for. 

The theoretical calculations, on the other hand, only 
contain the electronic contribution. Since the broaden- 
ing due to energy resolution is already included in the 
fitting procedure, the other experimental corrections are 
due to electron-phonon broadening (r e _ p /j) and final- 
state broadening (Tf), where r=2ImX! is the scattering 
rate. In general, all these contributions add linearly and 
give a total scattering rate F t such that 



r t = r c 



r, 



-ph 



T 



imp 



f 



(5) 



where r e _ e indicates the contribution to the linewidth 
due to electron-electron scattering. In the Debye model 
and for the high temperature limi t 53 ' , the electron- 
phonon broadening depends linearly on the temperature 
T as: 



e—ph 
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2n\k B T 



(6) 



where ks and A are the Boltzmann and electron- 
phonon coupling constants, respectively. Eq. ([6]) holds 
for temperatures T > Toebye/^ (here T=300 K and 
Ir>ebye(Co)=385 K). From resistivity measurements^ we 
obtain for Co A « 0.3 and as a result lmE e - p h « 25 meV. 



The final state broadening was determined by unfold- 
ing the bands around the A point and calculating the 
experimental electron initial and final state velocities as- 
suming a free electron parabola as final state. Firstly, the 
initial state velocities were calculated from the derivative 
of the E(kj_) experimental dispersions of the different 
majority and minority spin bands. Secondly, the inelas- 
tic mean free path values were taken from the semiem- 
pirical estimation given by Penn^i and available refer- 
ence datai£2 In general, due to the relatively flat dis- 
persions of the d bands, a maximum ratio between the 
initial and final state group velocities \vi±/vfj_\ ~ 0.04 
was obtained, equivalent to a maximum contribution to 
lmE exp of about 40 meV. This indicates that final-state 
effects are not the main broadening mechanism, as men- 
tioned above. Finally, the electron-impurity scattering, 
an energy and temperature-independent quantity, can be 
directly estimated as an average for both spin channels. 
It corresponds to the value of ImS eip exactly at the 
Fermi level position, and by extrapolation it amounts to 
luiEi mp ^50 meV. 

Although no quantitative agreement between the ex- 
perimental and theoretical ImE is obtained for Co, both 
spin channels are qualitatively reproduced in the theo- 
retical results. Interestingly, the experimental data sup- 
port the existence of a pronounced spin-dependent effect 
in the quasiparticle lifetimes in agreement with previous 
observations^ However, this is the first time that the use 
of spin resolution has been exploited to demonstrate such 
a notorious effect. In the experimental results, we obtain 
to a first approximation a k-independent IraT,Exp which 
exhibits an almost linear energy- dependent behaviour ap- 
proximated as Im£jj XJ) ~ 0.84E and ImK^ p ~ 0.16E 
for majority and minority spin electrons, respectively. 
Furthermore, the relative overall agreement between the 
experiment and the theoretical results indicates that in 
the case of Co, the role of non-local correlations is not 
as important as for Fe^ This could be attributed to a 
reduction of the non-locality of correlation effects with 
increasing atomic number. 

Regarding the theoretical calculations of ImS u mft in 
Fig. 7, it can be seen that increasing U does not strongly 
affect the global behaviour of the calculated curves. Note 
that this is something expected following our previous 
discussion concerning the shifting of the peaks with in- 
creasing U, which was also not very pronounced. In this 
respect, it should be emphasized that our results do not 
imply a violation of the Kramers-Kronig consistency^ 3 - 
which holds between the real and imaginary parts of the 
self-energy, neither in the experiment nor in the theo- 
retical calculations. This is because the Kramers-Kronig 
consistency results from a global transformation over an 
extremely large BE range, which runs up to infinity. This 
means that in the narrow BE energy interval studied 
here, the deviations between the experimental and theo- 
retical self-energies may become substantial without any 
noticeable effect in the global behaviour of the Kramers- 
Kronig transformation. To check this more in detail ex- 



12 



perimentally, time-consuming measurements over an ex- 
tremely wide BE range would be necessary. In Fig. 7, 
the deviations between experimental and theoretical data 
seem to be equally pronounced for both majority and mi- 
nority spin electrons, with the exception of the minority 
spin ImE dm ft for U=3 eV, which seems to agree better 
if the electron-impurity scattering contribution is sub- 
tracted from the experimental curve. Nevertheless the 
overall agreement shows that the experimental ImS is 
roughly a factor of 2 too large in average when compared 
to the calculations. In the LSDA+DMFT calculations, 
this discrepancy is more pronounced in the majority spin 
channel at higher binding energies than at lower ones. 
Note that both the IuiEdmft and I1T1S3B5 agree well 
for the same U values, with the exception that at higher 
binding energies they start to deviate from each other. 
We attribute this difference to the fact that LSDA+3BS 
calculations explicitly include electron-hole pair excita- 
tions which increase ImS. The extent of the present 
agreement contrasts a very recent non spin-resolved anal- 
ysis from Moulazzi et ali^ 4 -, in which excellent agreement 
between experiment and previously reported results of 
ImYiDMFT^ is achieved for both spin channels, most 
probably due to an arbitrary linewidth assignment of 
nearly-overlapping spin contributions. 

V. CONCLUSION 

In summary, we have shown a detailed comparison be- 
tween spin and angle-resolved photoemission experiments 
and state-of-the-art theoretical calculations of ferromag- 
netic hep Co(0001). We have analyzed in detail the agree- 
ment between theoretical many-body 3BS and DMFT 
based calculations and experimental spin-resolved pho- 
toemission spectra. The theoretical methods lead to 
an important improvement compared to LSDA calcula- 
tions by including, in an unified picture, many-body cor- 
rections with multiple-scattering, matrix-elements and 
surface-related effects. However, we have demonstrated 



that these theories, at least in the current implemen- 
tation, do not find exact quantitative agreement with 
the experiment, in particular concerning the scattering 
rates. Although further theoretical work is needed to 
determine the origin of these differences, a possible ex- 
planation could be related to the existence of non-local 
correlation effects. Such non-locality is caused by long- 
range electron-electron interactions, which are excluded 
from the Hubbard model. These non-local correlation 
effects would contribute with an additional k-dependent 
correction to the self-energy functions Sdmft or E3BS. 
This would result in second-order variations in the cal- 
culated ReS and ImE, which may lead to a much bet- 
ter quantitative agreement between theory and exper- 
iment. The fact that non-local correlation effects can 
have a certain influence on the intermediate energy states 
is closely related to what we found for Fe.— However, 
the better overall agreement between theory and experi- 
ment obtained for Co at intermediate U values indicates 
that non-local correlations become weaker with increas- 
ing atomic number. It would be very interesting to ap- 
ply similar investigations to itinerant magnets like the 
series FeSi, MnSi and CoSi, where we also expect the 
non-local character of correlations effects to play an im- 
portant role. Although non-local interactions are beyond 
the many-body theories used in this work, a promising 
method to overcome this limitation in the near future is 
the dual fermion approach^ Alternatively, it would be 
an important step to implement into realistic electronic 
structure calculations the recently proposed parameter- 
free extended DMFT+GW scheme 64 , in which both the 
on-site and off-site correlations are included. 
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